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We investigate dynamical tunneling in many dimensional systems using a quasi-periodically mod- 
ulated kicked rotor, and find that the tunneling rate from the torus to the chaotic region is drastically 
enhanced when the chaotic states become delocalized as a result of the Anderson transition. This 
result strongly suggests that amphibious states, which were discovered for a one-dimensional kicked 
rotor with transporting islands [L. Hufnagel et a/., Phys. Rev. Lett. 89, 154101 (2002)], quite 
commonly appear in many dimensional systems. 
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According to the semiclassical eigenf unction hypoth- 
esis, each eigenfunction of a quantum system, whose 
classical counterpart exhibits mixed- type phase space, 
is expected to be localized on either regular or chaotic 
components of classical phase space in the semiclassical 
limit [1 J. In a finite h regime, however, the tunneling 
effect comes in and it hybridizes localized eigenfunctions 
together. This has invoked considerable interest in quan- 
tum tunneling in two or more degrees of freedom sys- 
tems [2H5]. 

Since the tunneling effect is supposed to be exponen- 
tially small in general, one may regard tunneling merely 
as a correction to the semiclassical eigenfunction hypoth- 
esis. However, it was found that under a certain condi- 
tion eigenfunctions are not necessarily localized on either 
regular or chaotic regions even when h is much smaller 
than the area of the regular region. Thus it looks to be 
violating the semiclassical eigenfunction hypothesis [7J|8]. 

Such states are called amphibious eigenstates and their 
origin was discussed in Ref. [8|: They pointed out the 
importance of the relation between two time scales, the 
Heisenberg time Th in the chaotic sea and the torus's life 
time Tl, to understand such exotic states [9 j. Since the 
former is associated with the localization length in the 
chaotic sea, the argument on the two competing time 
scales can be rephrased using the interplay between dy- 
namical localization [10] and dynamical tunneling. 

A more direct explanation for the "flooding" of the 
wavefunction [TJ [8] can be given based on the time- 
domain semiclassical picture [11]. In particular, the the- 
ory of complex dynamical systems tells us that there are 
exponentially many complex orbits which connect the 
torus and chaotic regions. It could further be proved 
that they (i) act as tunneling orbits when they stay in 
the torus region, and (ii) behave as if they are real or- 
bits after reaching the chaotic sea [11 . In other words, 
dynamical tunneling and dynamical localization are con- 
trolled by the same orbits with amphibious characters, 
i.e., both regular and chaotic ones. Therefore, destruc- 
tive interference among the semiclassical orbits causing 
the suppression of classical diffusion [12] simultaneously 



inhibits chaotic tunneling from the torus to the chaotic 
sea. On the contrary, the attenuation of dynamical local- 
ization and the recovery of the chaotic tunneling occur 
simultaneously [15] . 

Amphibious states have been considered to appear in 
quite a specific situation where the system has acceler- 
ator modes [7] [8]. However, the accelerator modes are 
introduced just to prepare extremely large localization 
lengths [8]. The system with large localization lengths 
can also be realized in many degrees of freedom systems 
according to the correspondence of the dynamical local- 
ization to the Anderson localization [T4HT7] . This natu- 
rally leads to a conjecture: The appearance of amphibious 
states is a common feature in many- degrees of freedom 
systems. This, however, needs to be scrutinized, since 
Backer et a/.'s argument in Ref. [8] provides only a suffi- 
cient condition of the absence of amphibious states under 
the small localization lengths, and no condition that en- 
sures the amphibious states is known. In this Letter, we 
will provide numerical evidence of the conjecture. 

We examine a one degree of freedom kicked rotor whose 
phase space (q,p) is divided into a torus region p < b and 
a chaotic sea p > 6, which are connected by dynamical 
tunneling. To introduce many-dimensionality in effect, 
the system is described by the kicked Hamiltonian with 
quasi-periodic modulations p~5j [19] : 



H = T{p) + V{q)Y J 5{t-n) 

n 
1 M 



j=l n 

where V(q) = /ccos(27rg)/(47r 2 ), and modulation frequen- 
cies {uj/(2tt)}J^ 1 are irrational and non-resonant with 
each other [20]. T(p) and g(p) are specified in the follow- 
ing. We impose periodic boundary conditions < q < 1 
and -Wp/2 <p< W p /2 [21]. 

Before studying the mixed case, we first examine the 
fully chaotic case to identify insulator and metallic phases 
in the parameter space of e, with a sufficiently large value 
of W p . We choose T(p) = s(p — b) 2 /2 and g(p) — p — b. 
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FIG. 1: (Color online) (a) The scaling parameter £, which 
is proportional to the localization length and the inverse of 
the diffusion constant in the insulator and metallic phases, 
respectively, as a function of the strength of the modulation 
e in our model [Eq. (lib] . T(p) and g(p) are chosen so as 
to make the whole phase space chaotic (see the main text). 
£ diverges at the critical point e c — 0.227. A thin curve 
indicates £ = a\e — e c \~ u with v — 1.55 and a = 2.4 x 10 -2 . 
The parameters are set to be k = 2, s = 4 + v^/10, 6 = 0, 
M = 2, cji = 2tt^5, cj 2 = 2ttV13, W p = 1000, and Planck's 
constant h — 0.2. (b) The decay rate 7 of P£ [Eq. ([3j] as 
a function of e. Different marks correspond to ten different 
values of s for a small interval 4.17 < s < 4.19. For e<e c , 
chaotic tunneling is almost forbidden, i.e., 7^0 [18]. Around 
e rsj e c , 7 exhibits strong fluctuation, which becomes smaller 
for larger e. The parameters for Eq. ^ are 6 = 5.6 — W p /2, 
6 = 6 + 0.2, cj = >/3 and (3 = 50. The initial state |^o) 
satisfies the EBK condition whose center of the momentum is 
at p = b — 0.6 (see, Appendix of Ref. [13]). 



In the absence of modulation, i.e, e = 0, the mapping 
derived by Eq. (fTl) can be reduced to the standard map- 
ping whose nonlinearity parameter is K = sk. In the 
following, we examine the case that the nonlinearity is 
sufficiently large K ~ 8, where the phase-space is mostly 
filled with a chaotic sea [22]. When e is non-zero, this 
model is essentially the same as Casati et a/.'s quasi- 
periodic kicked rotor, which is equivalent with a M + 1 
dimensional tight-binding model with pseudo-disorder 
and is shown to exhibit the Anderson transition with 
M = 2 [151 [T7]. We also examine the case M = 2 and 
numerically confirmed that our model is in insulating 
and metallic phases at e < e c and e > e c , respectively, 




FIG. 2: (Color online) Poincare section around the torus re- 
gion of the tunneling degrees of freedom under the absence 
of modulation (i.e., e = 0). The torus region p < b and the 
chaotic sea p > b are seen to be sharply divided. The density 
plot of the Husimi function of the initial torus state, which 
is constructed by EBK quantization, is superposed. The pa- 
rameters are same as in Fig. Il[b), except s — 4 + v^/10. 



where e c ~ 0.227 is the metal-insulator transition point. 
Using the finite size scaling technique [TTJ [23] [24], we 
obtained numerical evidence of the Anderson transition 
in the chaotic sea. The details will be reported else- 
where [25]. We show how the scaling parameter £ [24] . 
which is proportional to the localization length and the 
inverse of diffusion constant in the insulating and metallic 
phase, respectively, depends on e in Fig. Ilia). 

Next, we introduce a torus region in p < b with keeping 
untouched the chaotic nature of the region p > b (see, 
Fig. [2]) , using the following kinetic term [13] 



T(p) = -(p-b) 2 e ( P -b)+co( P -b), 



(2) 



where 9p(x) — [1 +tanh(/te)]/2 is a smoothed step func- 
tion with a smoothing parameter (3 [26 . The region 
p < b is filled with tori because T(jp) is effectively linear 
there. At the same time, we employ the modulation term 
g{jp) = (p — b)9p(p — 6), where b is slightly larger than b. 
This makes the torus region almost independent of the 
modulation. Hence, only dynamical tunneling induces 
transitions between the torus region and the chaotic sea. 
Tunneling leakage from the torus to chaotic regions is 
monitored by the integration of |(p|^n)| 2 f° r the whole 
torus region [27] 



W p /2 



\(p\^n}\ 2 d P , 



(3) 



where \i/j n ) denotes the state vector at time step n. We 
prepare the system initially to be in a torus state |^o), 
which satisfies the Einstein-Brillouin-Keller (EBK) quan- 
tization in the torus region (see, Fig. [2]). Time evolution 
of P^ strongly depends on e, as shown in Fig. pi When 
e is far below e c , P^ keeps almost unity within, say, 10 6 
steps. As shown in Fig. J4^a), the corresponding wave 
function in the momentum representation exhibits dy- 
namical localization in the chaotic sea. As was shown 
in [28], and will be reported in [25], there appear expo- 
nentially many complex orbits connecting the torus and 
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FIG. 3: (Color online) The survival probability P^ of a torus 
state as a function of the number of kicks n. The parameters 
are the same as in Fig. [2] When the chaotic sea is in the insu- 
lating phase e < e c , the tunneling decay is almost forbidden. 
For e > 0.3 (> e c ), P^ exhibits exponential decay. Around 
the transition point, the decay of P„ strongly depends on e, 
and the decay is too slow to determine whether or not P„ 
obeys the exponential law. 



chaotic region, nevertheless such "flooding" of complex 
orbits are suppressed because of the destructive interfer- 
ence in the chaotic region. 

For larger e, chaotic tunneling recovers (see, Fig. [3]). 
Figure [lib) shows the e-dependence of decay rate 7, 
which is obtained by fitting of P^ . Around e c , 7 changes 
significantly with large fluctuations. A strong correlation 
between e dependences of 7 and £ is evident. Indeed, 
when 7 is significantly different from zero, the classical 
diffusion in the chaotic sea p > b is recovered, as shown 
in Fig. pFb). This strongly suggests that chaotic tunnel- 
ing recovers when the chaotic sea is in the metallic phase. 
Note that the localized component on the torus region al- 
most disappears at n = 10 6 , and the tails of wavefunction 
is clearly given as Gaussian (see inset of Fig. Kb)). 

We remark on the fluctuation of 7 (see, Fig. [ljb)). 
Around e c , 7 strongly depends on s as well as e. A pos- 
sible origin of the fluctuation is the resonance induced 
by near degeneracies among approximate quasienergies 
of torus and chaotic states [13]. Such a strong fluctu- 
ation against parameter variations is a common feature 
of dynamical tunneling in nonintegrable systems [2 . For 
smaller e (<C e c ), the resonances are ineffective [29] . In 
contrast, the effect of the resonances becomes prominent 
around e c . This is because the effective density of states 
around the torus state become larger due to the expo- 
nential growth of the localization length in the chaotic 
sea (see, Fig. []Ja)). For much larger e (^> e c ), the fluc- 
tuation becomes smaller. This is supposed to be due to 
the completion of the transition to the metallic phase. 
In other words, the effective density of states around the 
torus state become so large, the fluctuations induced by 
the resonances are averaged out. 

So far, we have focused on the case that W p is suffi- 
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FIG. 4: (Color online) Snapshots of the momentum distribu- 
tion \{p\i/j n )\ 2 . (a) In the insulating phase (e = 0.2), Kpl^n)! 2 
remains almost unchanged in the torus region and dynami- 
cal localization occurs in the chaotic sea. (b) In the metallic 
phase (e = 0.8), \{p\ipn)\ 2 decays in the torus region and the 
envelope of K^l^n)! 2 in the chaotic sea around the torus re- 
gion is almost unchanged. However, as is seen in the inset 
of (b), the diffusion toward the uniform distribution in the 
momentum space is not completed for n < 10 6 . 



ciently large. As for the "phase transition" between the 
suppression and restoration of chaotic tunneling, our nu- 
merical result suggests that the convergence for the "ther- 
modynamic" limit W p — » 00 is quite fast. In this ideal- 
ized limit, one may regard that the delocalized chaotic 
sea plays a role of "particle bath" . 

To be precise, however, the transition from the sup- 
pression of chaotic tunneling to the recovery is not sharp, 
but rather smooth, as seen in Fig. [ljb). This is because 
the recovery occurs even in the insulating phase with 
sufficiently large localization length due to the Backer et 
a/.'s condition Tl — Th [8 , which determines the border 
between the suppression and the recovery. 

Furthermore, the crossover between the suppression 
and the restoration of chaotic tunneling occurs even when 
W p is rather small. If e is sufficiently larger than e c , the 
time evolution of P^ obeys the irreversible decay for a 
short time period, and after that, exhibits erratic oscilla- 
tions (see, Fig. I5J. Although the physical picture based 
on the thermodynamic limit is inapplicable anymore, our 
numerical result indicates the presence of a considerable 
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FIG. 5: (Color online) Evolution of P^ for the case of smaller 
W p = 30. We choose e = 0.5 (> e c ). To ensure the spectrum 
of the system in the de localized "phase" to be discrete one, 
the modulation frequencies are chosen to be rational [15], i.e., 
wi/(27r) = 682/305 and cj 2 /(2tt) = 11/3 (x), 119/33 (□), 
649/180 (O), and 4287/1189 (A). Other parameters are the 
same as in Fig [3] 



number of amphibious states, i.e., eigenstates that have 
significant overlap with both torus and chaotic regions. 

We have numerically investigated the stability of torus 
states surrounded by chaotic seas in many degrees of free- 
dom systems. If the coupling strength exceeds a criti- 
cal value at which the Anderson transition occurs in the 
chaotic region, the nature of dynamical tunneling drasti- 
cally changes. The result of the smaller W p case may give 
rise to a reexamination of quantization condition of many 
dimensional mixed systems. Although the quantization 
of tori and chaotic seas can be carried out separately in 
the semiclassical limit, this seems not to be the case even 
when the size of Planck's constant is considerably smaller 
than the area of a torus region, due to the emergence of 
amphibious states. The "nonseparability" among regu- 
lar and chaotic region certainly has been a problem of 
tunneling in nonintegrable systems. Our result suggests 
that the occurrence of the nonseparability, i.e., the emer- 
gence of amphibious states, is the rule rather than the 
exception in many-dimensional mixed systems. 
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